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^ A geometric  optics  analysis  is  made  for  a wavefront  sensor  that  uses  a 
moving  Ronchi  grid.  It  is  shown  that  by  simple  data  processing  a signal 
proportional  to  the  gradient  of  the  distorted  wavefront  can  be  obtained  for 
each  element  of  an  array  of  photodetectors  that  cover  the  entire  image  oi  a 
telescope  aperture.  The  effect  of  diffraction  is  calculated  for  this  system. 
This  shows  the  limitations  and  provides  a justification  for  the  geometric 
analysis, 
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I.  INTRODUCTION 


There  are  several  types  of  adaptive  optical  systems  being  considered 
1-3 

or  being  developed  for  imaging  an  object  through  a turbulent  atmosphere. 
Some  of  these  use  a wavefront  sensor  to  measure  the  degree  of  wavefront 
deformation  caused  by  atmospheric  turbulence.  This  information  is  then  used 
to  deform  the  surface  of  a flexible  mirror  or  control  some  other  active  ele- 
ment to  compensate  for  the  wavefront  deformation  and  thereby  improve  the 
image  quality.  The  wavefront  sensor  used  in  one  of  these  adaptive  optical 
imaging  systems,  which  has  been  described  by  Feinleib  and  Hardy,  ^ consists 
of  a rotating  Ronchi  grid  placed  in  the  image  plane  of  the  telescope.  This 
wavefront  sensor  is  relatively  simple  and  insensitive  to  external  disturbances. 
The  operation  is  explained  by  Feinleib  and  Hardy  by  means  of  a physical 
optics  theory,  in  which  the  sensor  is  assumed  to  behave  as  a shearing  inter- 
ferometer. This  report  presents  a geometric  optics  analysis  of  this  wave- 
fror.c  sensor.  The  justification  for  the  geometric  optics  approach  is  given  in 
the  Appendix.  One  result  of  this  analysis  is  a simplification  in  the  data  pro- 
cessing required  to  translate  sensor  data  into  wavefront  deformation 
information. 


^Julies  Feinleib  and  J.  W.  Hardy,  "Wideband  Adaptive  Optics  for  Imaging,  " 
Proc.  S'PIE  75,  103  (1976). 

^A.  Buffington,  F.  S.  Crawford,  R.  A.  Muller,  A.  J.  Schwemin,  and 
R.  J.  Smits,  "Active  Image  Restoration  with  a Flexible  Mirror,  " 

Proc.  SPIE  75,  89  (1976). 

3 

Virendra  N.  Mahajan,  "Optical  Wavefront  Correction  in  Real  Time,  " 

Proc.  SPIE  75,  109  (1976). 
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II.  BASIC  THEORY 


Conaider  a simplified  model  of  a telescope  with  a moving  grid  wavefront 
sensor,  as  shown  schematically  in  Fig.  1.  A flexible  corrector  mirror  and 
an  image  sensor  normally  included  in  a compensated  imaging  system  are  not 
pertinent  to  this  discussion  and  are  not  shown  in  Fig.  1.  The  telescope  focuses 
a distant  object,  defined  by  the  angular  coordinates  (>,  P),  onto  an  image  plane 
(x,  y).  The  aperture  of  the  telescope  lies  in  a (u,  v)  plane,  and  the  deviation 
of  the  incoming  wavefront  from  a plane  wavefront  at  this  aperture  is  described 
by  A(u,  v).  It  is  assumed  that  the  angular  extent  of  the  object  is  s\ifficiently 
small  that  all  rays  from  the  object  to  a given  point  on  the  aperture  plane  ex- 
perience the  same  deviation  in  optical  path  length,  i.  e. , the  object  is  smaller 
than  the  aplanatic  patch  size.  For  a plane  wavefront,  the  image  formed  by  the 
telescope,  I(x,  y),  is  a true  mapping  of  the  object,  i.e.. 


Ij(x,y)  = Ijj(x/F,y/F) 


if  the  telescope  is  assumed  to  be  perfect. 

It  che  wavefront  is  distorted,  the  rays  of  light,  which  are  perpendicular 
to  the  wavefront,  are  deviated  by  small  angles  of  [8A(u,  v)/8u]  and  [8A{u,  v)/8v] 
m the  X and  y directions,  respectively.  Divide  the  aperture  into  a number  of 
equal  square  segments  of  area  ^A^^  centered  on  the  coordinates  m and  v.. 

Also,  make  the  elements  sufficiently  small  that  the  gradient  of  the  wavefront 
is  essentially  constant  over  each  element.  Each  of  these  elements  of  the 
telescope  aperture  then  forms  a complete  image  that  is  shifted  by[8A(u,  v)/ 
8u]_F  in  the  x direction  and  by  [8A(u,  v)/8  v], ^F  in  the  y direction.  The  ele- 
mental image  formed  by  each  AA. . is 


AIj-’Cx,  y) 


(1) 


PLANE 


Figure  1.  Telescope  and  Wavefront  Sensor 
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and  the  complete  image  is  a superposition  o£  all  ^l^^#  i*  e*. 


Ii(x.y) 

i j 


This  is  a noncoherent  point  of  view  for  which  it  is  assumed  that  the  intensities 
of  the  elemental  images  add  together  without  interference.  This  sxunmation, 
however,  could  just  as  well  be  a sum  of  the  complex  amplitudes.  The  signals 
on  the  photodetectors  or  image  sensors  are  then  proportional  to  the  '.bsolute 
value  of  the  square  of  the  amplitudes.  The  gradients  of  the  wavefront  dis- 
tortion in  the  x direction  are  determined  by  the  use  of  a Ronchi  grid,  with 
bars  oriented  parallel  to  the  y axis,  moving  in  the  x direction.  The  gradients 
of  the  distorted  wavefront  in  the  y direction  are  measured  by  a second  moving 
Ronchi  grid  oriented  at  a right  angle  to  the  first.  The  second  Ronchi  grid  is 
not  shown  in  Fig.  1.  The  light  transmitted  through  the  grid  is  collected  by  a 
lens  that  projects  an  image  of  the  telescope  aperture  onto  a matrix  of  photo- 
tube detectors. 

The  operation  of  the  wavefront  sensor  can  be  explained  by  considering 
only  one  grid,  which  can  be  defined  by  an  aperture  function  h(x),  such  that 


h(x)  = 1 for  -d/4  S x < d/4 
h(x)  = 0 for  -d/2  <x<-d/4  and  d/4  <x<d/2 


(2) 


where  h(x)  is  an  even  periodic  function  with  a period  d.  The  grid  is  translated 
with  a velocity  V in  the  x direction  such  that  the  transmission  of  the  grid 
becomes  h(x  - Vt).  Assume  that  d is  made  larger  than  the  image  size  so 
the  image  is  transmitted  through  only  one  slit  at  a time.  The  moving  grid 
function  is  shown  in  Fig.  2. 

The  matrix  of  photodetectors  corresponds  to  the  square  elements  of 
the  telescope  aperture  such  that  each  photodetector  measures  the  light 
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transmitted  through  a given  element  of  the  aperture.  The  light  incident 

on  a given  photodetector  is  given  by 


y)h(x  - Vt)  dx  dy 


(3) 


where  the  integral  is  taken  over  the  image  area.  With  the  use  of  Eq.  (1), 
we  obtain 


1.  i (^»  v)\ 


h(x  - Vt)  dx  dy 


(4) 


The  problem  is  to  extract  the  value  of  the  wavefront  gradient  with  a 
simple  data  processing  scheme.  Since  h(x)  was  chosen  to  be  an  even  function, 
we  multiply  w^^(t)  by  an  odd  function  and  integrate  over  one  period.  Define 
the  function  g(t)  to  be  a periodic  square  wave,  i.  e.  , let 


g(t)  = +1  for  0 i t < d/2v 
g(t)  = - 1 for  -d/2v  5 t < 0 


The  integral  over  one  period  is 


(t)  g(t)dt 


(5) 


This  integral  is  directly  proportional  to  the  gradient  of  the  distorted  wave- 
front  with  some  restrictions. 
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A simple  way  to  demonstrate  the  above  assertion  and  to  determine  any 
restrictions  is  to  consider  the  case  of  a square  image  of  size  a,  such  that 
a<  d/2. 


Ij(x)  = Iq  for  -a/2  £ x < a/2 

= 0 for  x<-a/2  and  x > a/2 

The  light  transmitting  through  the  grid  is  the  convolution  of  the  image 
distribution  and  the  traveling  grid  function.  From  Fig.  2 it  can  be  seen 
that  the  transmitted  light  is 


(t)  = Iq 


h(x  - Vt)dx 


(6) 


The  result  of  the  integration  for  this  case  is  a train  of  trapezoidal  pulses  as 
shown  in  Fig.  2.  This  temporal  function  is  then  multiplied  by  g(t)  and 
integrated  over  one  period 


^x=  f <*)g<t)dt 

T 


This  is  done  graphically  in  Fig.  2,  and  the  result  is 


(7) 
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Therefore,  the  integral  over  one  period  is  indeed  proportional  to  the  gradient 
of  the  distorted  wavefront.  If  the  gradient  is  zero,  the  integral  over  one 
period  is  zero.  In  order  to  extract  the  wavefront  gradient,  it  is  necessary  to 
divide  by  the  integral  of  the  transmitted  light  over  one  period,  namely 


4' 


which,  for  this  case,  is 


aa. 


-iJ.  4- 


"o  A 2V 


Then 


/5A(u,  v)\  _ 

\ 6u  “ * 4F 


d X 


4F 


(8) 


The  analysis  becomes  more  complex  if  the  flat  top  of  the  trapezoid  does  not 
extend  beyond  t > 0,  ie,  if 


tj 


In  the  latter  case,  quadratic  terms  are  introduced,  and  the  sensitivity  is 
decreased. 


Extension  of  the  above  analysis  to  an  object  of  arbitrary  shape  can  be 
made  graphically  with  the  information  given  in  Fig.  3.  The  preceding 
analysis  still  holds  if  the  transmitted  pvilse  is  approximately  symmetric, 
and  if  the  shifted  image  remains  within  the  size  of  the  grid  opening  at  t = 0, 
viz,  -d/4  < x < d/4. 


III.  SUMMARY  AND  CONCLUSIONS 


It  has  been  shown  that  the  operation  of  a moving  grid  wavefront  sensor 
can  be  explained  by  a simple  geometric  optics  analysis.  The  gradient  of  the 
distorted  wavefront  at  any  given  aperture  location  can  be  extracted  by  multi> 
plying  the  transmitted  light  by  an  odd  square  wave  function  and  integrating 
over  one  period.  A square  object  was  used  as  an  example  to  demonstrate 
the  operation  and  to  determine  the  approximate  limits  of  operation. 
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APPENDIX.  EFFECTS  OF  DIFFRACTION  IN  THE  WAVEFRONT  SENSOR 

The  validity  of  the  geometric  analysis  of  the  moving  grid  wavefront 
sensor  for  any  specific  case  requires  that  the  effects  of  diffraction  on  the 
operation  of  the  sensor  be  small.  The  most  important  effect  of  diffraction 
is  the  spreading  of  light  from  one  element  of  the  aperture  onto  photodetectors 
associated  with  adjacent  elements.  This  diffraction  effect  is  caused  by  the 
presence  of  the  Ronchi  grid  in  the  image  plane.  It  is  therefore  necessary  to 
calculate  the  fraction  of  light  diffracted  outside  the  geometric  image  of  the 
aperture  element  and  to  determine  the  spatial  period  of  the  Ronchi  grid  for 
which  this  diffraction  effect  is  acceptably  small. 

For  a telescope  of  diameter  D and  focal  length  F the  diameter  of  the 
point  spread  function,  the  Airy  disc,  at  the  focal  plane  is 


2.44X 

D 


F = 1. 342  X 10 


6 F 
D 


for  a wavelength  of  0.  55  ^im.  Now,  mask  off  the  telescope  aperture  except 
for  one  square  element  of  size  a,  where  a is  an  integral  fraction  of  the 
telescope  diameter.  The  amplitude  point  spread  function  from  this  element 
is  a sine  pattern  in  both  the  x and  y directions 


AU(x,  y)  = a sine 


where  k = 2it/X. 

The  size  of  the  central  lobe  of  this  diffraction  pattern  is  (2\/a)F.  If  N 
is  the  number  of  elements  across  the  telescope  diameter,  then  the  size  of  this 
pattern  is  approximately  N times  the  size  of  the  Airy  disc.  If  a screen  or 
image  sensor  were  placed  at  the  focal  plane  of  the  telescope,  all  of  the  images 
from  the  individual  elements  would  add  together  coherently  to  produce  the 
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Airy  disc;  but  in  the  absence  of  a recording  medium,  all  of  these  images  can 
be  considered  independently  because  of  superposition. 

Figure  1 shows  how  the  aperture  elements  are  imaged  onto  a phototube 
array  by  a lens  placed  behind  the  image  plane  of  the  telescope.  In  the  absence 
of  the  Ronchi  grid,  the  image  of  the  telescope  aperture  formed  by  this  lens 
would  be  quite  sharp  and  would  be  limited  only  by  the  diffraction  limit  of  this 
lens.  The  Ronchi  grid,  however,  is  a restrictive  aperture  that  causes  signi> 
ficant  diffraction  and  spreading  of  the  telescope  aperture  image.  The  geo- 
metric analysis  determined  that  the  aperture  width  of  the  Ronchi  grid  should 
be  larger  than  the  image  size.  A reasonable  object  size  is  10  (xrad,  which 

is  smaller  than  the  usual  size  of  the  aplanatic  patch.  The  period  of  the  Ronchi 

-5  -5 

grid  should  then  be  d > 2 x 10  F.  Let  d = 2.  5 x lO  F.  The  image  of  the 
aperture  element  formed  on  the  phototube  array  is  approximately  the  fourier 
transform  of  the  product  of  the  Ronchi  grid  function  and  the  amplitude  point 
spread  function  of  the  square  aperture  element.  Define  the  grid  aperture  by 

g(x,  y)  = 1 from  -d/4  x d/4  and 
g(x,  y)  = 0 for  X < -d/4  and  for  x < d/4 


A single  aperture  has  been  assumed  instead  of  a periodic  function  for  g(x,  y); 
this  can  be  accomplished  by  placing  a slit  aperture  of  width  d in  front  of  the 
Ronchi  grid.  We  can  therefore  avoid  the  complex  effects  of  partially  coherent 
light  transmitted  through  adjacent  grid  openings.  The  complex  amplitude  of 
the  light  just  after  transmission  through  the  Ronchi  grid  is  then  AU(x,  y)g(x,  y) 

= AU'(x,  y).  The  image  of  the  aperture  on  the  plane  of  the  phototype  array, 
which  is  designated  the  (S,T))  plane,  is  the  fourier  transform  of  AU*(x,  y).  This 
fourier  transform  is  the  convolution  of  the  fourier  transform  of  AU(x,  y), 
which  is  just  tlie  geometric  image  of  the  square  aperture  element,  and  the 
fourier  transform  of  g(x,  y).  The  latter  is  given  by 


G(?,T1) 


^d/4  5) 

•'-d/4 


1 


3 
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The  size  of  the  central  band  (between  the  first  zeroes)  of  this  sine  function 
is  4Xf/d  = 0.088  f/F  (meters),  and  the  size  of  the  geometric  image  is  af/F 
(meters).  A reasonable  size  for  an  aperture  element  is  the  coherence  width 
of  the  light  transmitted  through  the  atmosphere,  which  varies  from  5 to  15  cm. 
We  therefore  choose  a value  of  0. 1 m for  the  size  of  an  aperture  element.  In 
order  to  graphically  represent  the  diffraction  spreading,  in  Fig.  4,  we  plot 
the  convolution  of  the  sinc^  function  given  above  with  a square  element  for 
two  cases.  The  first  case  is  for  equal  widths,  which  is  approximately  the 
preceding  case,  and  the  second  is  for  an  element  image  size  twice  that  of  the 
central  lobe  of  the  sinc^  function,  which  is  equivalent  to  doubling  the  period 
of  the  Ronchi  grid.  A graphical  integration  shows  that  approximat''ly  20%  of 
the  light  is  diffracted  to  adjacent  elements  for  the  first  case,  and  12%  of  the 
light  is  diffracted  out  of  the  geometric  image  in  the  second  case.  An  analysis 
of  the  effect  of  this  diffracted  light  on  system  operation  should  be  made  to 
determine  the  acceptable  level  of  light  spilled  into  adjacent  elements.  It  is 
clear,  however,  that  a wider  grid  spacing  is  preferable.  The  diffraction  in 
the  first  case  may  well  be  marginal  or  even  unacceptable. 


(a) 


(b) 


Figure  4.  Diffraction  Spreading  of  Image  of  Aperture  Element 


